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ABSTRACT 


The  numerical  ice  prediction  system  described 
in  rp  8 has  been  in  use  since  1971.  Its  outputs, 
however,  have  been  related  only  to  ice  growth,  not 
ice  drift.  An  addition  to  the  system  now  provides 
an  output  of  daily  ice  drift  at  62  arctic  stations. 
Drift  is  calculated  as  a function  of  wind  stress, 
water  stress,  and  Coriolis  force.  The  computer 
output  is  sufficiently  accurate  for  ice  drift  in 
the  locality  of  the  station  whose  wind  observations 
were  input  to  the  model.  However,  more  recent 
research  indicates  that  for  greater  accuracy  and 
for  confutations  of  divergence  and  convergence 
interaction  between  floes  must  be  considered. 
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FOREWORD 


Fleet  operations  in  the  Arctic  are  often  hindered, 
and  sometimes  aborted,  by  the  severity  of  ice  conditions. 
The  Naval  Oceanographic  Office  has  pioneered  in  the 
development  of  ice  forecasting  support  designed  to  per- 
mit the  Navy  to  operate  effectively  at  all  times  in  arctic 
waters.  This  publication  describes  an  automated  program 
for  computing  ice  drift  for  selected  locations  in  the 
Arctic.  It  supplements  a previously  existing  program  for 
computing  and  predicting  ice  thickness  for  these  locations. 
Ice  drift,  as  well  as  ice  thickness,  information  is  re- 
quired for  mine  warfare,  submarine  surfacings,  surface 
ship  transits,  and  sonobuoy  plants. 


CONTENTS 

[. 

INTRODUCTION 

1 

El. 

DATA 

1 

III. 

METHOD 

1 

A.  Mean  Wind  Speed  and  Direction 

1 

B.  Drift  Computation 

3 

IV. 

VARIATION  FOR  CONCENTRATION  AND  ROUGHNESS 

9 

V. 

CONCLUSIONS 

9 

VI. 

DISCUSSION 

9 

REFERENCES 

12 

FIGURES 

1. 

Sample  output  of  wind  velocity  and  ice  drift 

2 

2. 

Flow  chart  of  ice  drift  subroutine 

4 

3a. 

Velocity  diagram 

16 

3b. 

Force  diagram 

16 

TABLES 

1. 

Factor  Table 

10 

APPENDIXES 

A. 

Subroutine  Listing 

13 

B. 

Theoretical  Development  of  Wind  Drift  Equations 

16 

, w.rj 

V 

c'C<, 

r *•>.>*! 

I .  INTRODUCTION 


Since  July  1971  a numerical  ice  prediction  system  has  been  operated 
by  the  Naval  Oceanographic  Office  to  compute  various  sea  ice  characteristics 
and  forecasting  aids  of  importance  to  arctic  operations.  Gerson  (1975) 
completely  describes  the  system  as  of  1975;  therefore,  it  will  only  be 
briefly  discussed  here. 

The  original  outputs  of  the  prediction  system  include  the  following: 

* present  degree-day  accumulations 

* forecast  degree-day  accumulations 

* daily  synoptic  temperatures 

* temperature  trends 

* selected  sea  surface  temperatures 

* initial  freezeup  forecasts 

* reported  snow  depths 

* ice  thickness  forecasts 

All  this  information  relates  to  ice  formation,  growth,  thickness,  and  decay 
in-situ  without  consideration  of  ice  drift.  Since  the  knowledge  of  ice 
drift  is  extremely  important  to  the  forecaster,  its  inclusion  as  part  of 
the  daily  output  was  considered  desirable.  The  purpose  of  this  report  is  to 
describe  this  significant  update  to  the  ice  prediction  system,  namely  the 
computation  of  wind  drift  of  ice  at  each  of  the  62  stations  in  the  present 
network.  An  example  of  a typical  output  listing  is  shown  in  figure  1. 

II.  DATA 


The  data  are  the  same  as  used  for  the  basic  system,  and  are  fully 
described  in  Gerson  (1975).  A slight  revision  in  the  heading  is  discussed 
in  U.  S.  Department  of  Commerce  (1973).  The  data  base  remains  the  entire 
world's  synoptic  weather  observations  collected  by  the  National  Oceanic  and 
Atmospheric  Administration,  but  observations  for  each  day  now  are  on  one 
tape  rather  than  two  and  are  processed  on  its  IBM  360/195  computer. 

III.  METHOD 


A.  Mean  Wind  Speed  and  Direction 

The  wind  drift  is  derived  entirely  from  values  of  synoptic  wind 
speed  and  direction  as  reported  by  the  network  stations.  The  four  (or  fewer 
if  reports  are  missing)  daily  synoptic  wind  observations  are  vectorially 
averaged  to  produce  a mean  wind  speed.  The  equations  used  are  as  follows: 


Figure  1. 

Sample  output  of  wind  velocity  and  Ice  drift 


I = i > VinDi 


l V*  v. 
k 2-f  1 


icosDi 


where 


D = arctan  (I/J) 


V = (i2  + J2) 


k = number  of  observations  (1-4)  during  the  day, 


= wind  direction  at  time  i. 


V.  = wind  speed  at  time  i, 


D = mean  wind  direction  for  the  day,  and 


V = mean  wind  speed  for  the  day. 


The  computer  program  accounts  for  the  quadrant  in  these  computations.  This 
subroutine  is  listed  in  appendix  A,  and  a flow  chart  is  shown  in  figure  2. 

B.  Drift  Computation 


The  method  used  for  ice  drift  computation  is  that  of  Shuleikin  (1953) . 
This  method  has  been  used  by  Wittmann  and  MacDowell  (1964)  to  produce  a set 
of  curves  for  manual  computation  of  ice  drift.  In  order  to  obtain  a daily 
output  for  62  stations,  however,  an  automated  method  is  necessary.  The 
equations  as  presented  by  Shuleikin  do  not  immediately  lend  themselves  to 
computer  programming,  and  considerable  transformation  was  required.  A complete 
derivation  of  the  equations  used  in  the  computer  program  is  given  in  appendix  B. 
The  equations  used  are  as  follows: 
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-*fp*«w**  f-  - 


SET  INITIAL  VALUE  FOR  Z 


l/COS(PSI  (I)  ) **  Z 


DETERMINE  FINAL  VALUE 
BY  ITERATION  ON  Z 


STOP 


Z(I  = I)  = SQRT(Z(I) ) / (XWMNS (J)  * 1.68  **  Z/C2  **  Z + 
2 * Cl  * Z (I)  **  . 2 5 /XWMNS ( J) /CZ/I .68  + Cl  **  Z + 1 


PSI (1+1)  = ACOS ( 1 /SQRT ( Z ( I + 1))) 


1DRD(J)  = IWMND(J)  1 PSI(I  + 1) 


Figure  2. 

Flow  chart  of  ice  drift  subroutine  (con.) 


psi  (i  + \)Vy 

YThs 

CALCULATE  TRUE  DRIFT  ANCLE 

PSI  (1)  IS  FINAL 
DRIFT  ANGLE  RELATIVE 
TO  WIND  DIRECTION 

CALCULATE  DRIFT  SPEED 


XDRS(J)  = XWHNS(J)  * SQRT ( . 004  * SQRT(Z)  * PI/ (473  * F))  * 
COS (27)  + COS (18)  * SQRT(. 000625  * COS(PSI(I  + l))/COS(18)) 


Figure  2 . 

Flow  chart  of  ice  drift  subroutine  (con.) 


(1)  drift  direction 


i p = arc  cos 

z - » ?£l  2^  + CjJ  + 1 

v2c2!  vc2 


f = 1.458  x 10  4sin(lat)  ; and 


where 

V = wind  speed  (ft/sec)  , 

\p  = drift  to  the  right  of  the  wind  direction 
u = drift  speed  (ft/sec) 
lat  = latitude , and 

I = ice  thickness  (inches) 

The  equation  for  Z is  required  to  be  in  an  iterative  form,  thus  it  is  used  as 
follows : 


IV. 


VARIATION  FOR  CONCENTRATION  AND  ROUGNHESS 


The  Shuleikin  equations  assume  a single  ice  floe  and  do  not  consider 
surface  roughness.  An  empirical  table  of  wind  factors  (Zubov,  1945)  was 
compiled  by  P.A.  Gordienko  from  many  observations  of  ice  drift  in  the  Soviet 
Arctic.  Comparison  of  this  table  with  the  output  of  the  Shuleikin  equations 
shows  a match  at  0.5  oktas  and  4 tenths  extent  of  ridging  and  hummocking. 

A factor  table  (table  1)  was  then  constructed  for  use  in  correcting  the 
computer  output  to  the  applicable  concentration  and  roughness.  Instructions 
for  using  the  table  are  given  beneath  the  table.  As  an  example,  suppose  we 
have  a mean  wind  speed  and  direction  of  9.2  kn  and  102  deg,  which  produce 
an  ice  drift  and  direction  of  .34  kn  and  123  deg  (see  station  70086  in 
figure  1).  If  the  concentration  were  0.5  oktas  and  the  extent  of  ridging  and 
hummocking  4 tenths,  then,  referring  to  table  1,  the  ice  drift  speed  would 
remain  .34  kn  since  the  multiplicative  factor  is  1.  If,  however,  the 
concentration  was  6 oktas  and  the  ridging  3 tenths,  the  ice  drift  speed  would 
be  .10  kn  since  the  multiplicative  factor  is  . 3 . In  either  case,  the  drift 
direction  remains  123  deg,  since  this  table  only  relates  to  drift  speed. 

V.  CONCLUSIONS 


The  methods  described  above  have  been  programmed  and  are  now  operational 
on  the  UNIVAC  1108  computer  at  the  Naval  Oceanographic  Office.  Ice  drift 
information  is  now  regularly  produced,  together  with  the  other  products 
described  in  the  Introduction.  The  outputs  have  been  evaluated  by  comparing 
them  first  with  the  manual  computation  method  of  Wittmann  and  MacDowell  (1964) 
and  second  with  manual  calculations  accomplished  with  the  original  equations. 
These  comparisons  have  proved  the  computer  program  to  be  reliable.  Evaluation 
with  actual  observed  drift  conditions  has  not  been  accomplished,  but  since  the 
curves  from  Wittmann  and  MacDowell  have  been  used  successfully  for  years, 
the  computer  output  is  assumed  accurate. 

VI.  DISCUSSION 


As  previously  stated,  this  computation  of  ice  drift  is  based  on  Shuleikin's 
work  of  1953.  This  model  was  chosen  since  it  had  been  used  manually  for  years 
and  had  been  found  to  give  accurate  results.  In  addition,  it  was  most  easily 
adapted  to  estimating  wind  drift  from  station  observations  that  were  available 
in  the  existing  system. 

A basic  problem  exists  with  the  Shuleikin  model  in  that  it  only  considers 
the  motion  of  a single  floe  as  it  results  from  wind  and  water  stress  and 
Coriolis  force.  No  consideration  is  given  to  the  effect  of  interaction 
between  floes.  In  the  program  herein  described,  the  effects  of  concentration 
and  ridging  are  considered  empirically  but  are  not  input  to  the  theoretical 
model.  Considerable  research  has  been  done  since  Shuleikin's  work.  Many  of 
the  recent  papers  are  published  in  the  extensive  AIDJEX  series  (University  of 
Washington,  1972-  ).  Campbell  and  Rasmussen  (1972)  begin  with  an  excellent 


-9- 


Table  1 


FACTOR  TABLE 

TOTAL  CONCENTRATION  OF  ICE  (OKTAS) 


0.5 

1 

2 

3 

4 

5 

6 

7.0 

7.5 

E* 

1 

.26 

.24 

.20 

.17 

.13 

.10 

.09 

.07 

.07 

2 

.51 

.47 

.40 

.34 

.28 

.22 

.18 

.15 

.12 

3 

.75 

.71 

.63 

.55 

.47 

.39 

.30 

.22 

.19 

4 

1.00 

.96 

.84 

.72 

.62 

.51 

.39 

.29 

.24 

5 

1.26 

1.19 

1.06 

.92 

.79 

.65 

.52 

.38 

.31 

6 

1.51 

1.43 

1.26 

1.10 

.94 

.80 

.63 

.46 

.38 

7 

1.75 

1.66 

1.47 

1.29 

1.10 

.92 

.73 

.55 

.46 

8 

2.01 

1.90 

1.69 

1.49 

1.28 

1.07 

.85 

.64 

.53 

9 

2.26 

2.13 

1.90 

1.64 

1.40 

1.17 

.94 

.71 

.60 

•Extent  of  ridging  and  hummocking  (tenths) 

Table  of  multiplicative  factors  for  estimating  speed  of  ice  drift 
from  computer  printout.  The  drift  speeds  given  in  the  printout  are 
true  for  0.5  oktas  and  4 tenths  ridging.  To  determine  correct  speed 
enter  concentration  and  extent  of  ridging  into  this  table  and  multiply 
computer's  drift  speed  by  the  factor  thus  obtained. 
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✓ 


summary  of  ice  dynamics  up  to  1972  and  proceed  to  give  a model  that  considers 
ice  as  a viscous  material.  Rothrock  (1975)  considers  sea  ice  as  a plastic 
material.  Additional  papers  are  expected  shortly. 

Ice  drift  estimates,  computed  by  the  model  discussed  in  sections  III  and 
IV  of  this  paper,  are  sufficiently  accurate  for  use  in  the  locality  of  the 
station  whose  wind  reports  are  used  in  the  c:  lculation.  They  are  not 
sufficiently  accurate  for  the  calculation  of  areas  of  convergence  and 
divergence  or  for  interpolation  between  stations  or  extrapolation  into  areas 
lacking  wind  observations.  To  add  these  features  to  the  system  in  the  future, 
it  will  be  necessary  to  implement  one  of  the  more  sophisticated  models, 
requiring  more  specific  inputs  not  now  readily  available. 
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APPENDIX  B.  THEORETICAL  DEVELOPMENT  OF  WIND  DRIFT  EQUATIONS 

The  frictional  stress  between  two  media  is  assumed  to  be  proportional 

to  the  square  of  the  relative  velocity  between  them.  Let  be  the  wind 

velocity,  V±  be  the  ice  velocity,  assumed  to  move  in  +x  direction  and 

V be  the  water  velocity.  The  frictional  stress  (T  ) between  air  and 
w a 

ice  is  then  given  by: 

|Tal  = ValV^il2  = kaPa|Va'2  since  ^a1  >>  K' 


The  stress  (T±)  between  ice  and  water  is  given  by: 


lTil  = WVi  'V 

In  addition  to  these  two  forces  we  have  the  Coriolis  force: 

-►  “►  “► 

C = f (Vixk) 

acting  on  the  ice  where  f = 2wsine.  The  motion  is  assumed  to  take  place  under 
the  equilibrium  of  these  three  forces  therefore  they  must  vectorially  add  to 
zero. 

We  now  choose  a coordinate  system  whose  X - axis  coincides  with  the 
direction  of  V , thus  producing  figure  3 . 


T tC-T^  =*  0 


Figure  3a  Velocity  diagram 


Figure  3b  Force  diagram 


i 


Referring  to  figure  3a  we  see 


, BEST  AVAILABLE  COPY 


IVl-vJ 


_ sin(ir/4-a) 
sin (a) 


_ sin  (tt/4)  cos  (a) -sin  (a)  cos  (tt/4) 

sin (a) 

= sin  (tt/4)  cot  (a) -cos  (tt/4)  = t, 

cot  (a) -cot  (tt/4)  = cos  (a) -1  = ^ 

sin  (tt/4) 

cot (a)  = 1+  2C 

again  from  figure  3a  we  have : 


V.  = V cos(a)+  V.  -v  cos(ir/4-a) 
l w 1 w 

Now  referring  to  figure  3b: 

T.  = T.*i  = k p |v.-V  I 2cos  (TT/4-Ot) 
ix  l w w 1 w 


'•  V! 


T.  = T . • j = k p V.-V  sin(TT/4-a) 
ly  l w w i w1 


T =T*i  = kp  Iv^|2cos4' 
ax  v a a a 


i ',y  ’’  wcw^ 


SC  aos+(c)  aoo?l 


bne 


I V-  . V I Q A 

1 W Z 1 W Vv 


T = T . j = k p V 2sin>F 
ay  v a a a 


: 9VBC  9W 


C = C*i  = -mf(kxV.)*i  = o=-mf (V.xk) *i 
x i 1 


C = C*J  =-mf(kxVi),j  = -mf|vi|  = f^xkj'j 


where:  m = mass  per  unit  area 


— + (»~fr\:r)  rcsd  = H'ncJ 


: sisriw 


*C  (C-fr\TT)  aoc+eaoor  n3m 


(b-i\Tr\  aocMQ.T>! 

k = the  coefficient  of  friction  between  ice  and  water  * ' 

w 


k = the  coefficient  of  friction  between  air  and  ice 
Pw  = the  density  of  water  $..v|  A - 


: sonic  s-oV 


, V-  -V  - .T  - son  T I T 

w j:  1 ■>.  v si  'v  1 xv 


p = the  density  of  air 


i i k = the  unit  vectors  along  the  x,  y,  and  z directions  respectively 

» j » » ■.  i i - 
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therefore: 


EY  forces  = -kwpj  Vj-vJ  2sin  (n/4-a)  -mf  | v±  | + |Tv|  sin'i'  = 0 


Ex  forces  *=  -kwPw|v^-Vw|  2cos (TT/4-a)  |tv|  cos’?  ■ 0 


Solving  for  sinH*  and  cosT  and  taking  the  ratio: 


tan4*  = 


I Vi  I +kwPwl Vj~Vwl  asin(TT/4-a) 
kwPwl^i"^wl  2c°s  (lr/4-a) 


= tan(TT/4-a)  + 


mf  | Vj  I 

\pJW\2cos^/4-cx) 


Now  since: 


IVi-vJ 


Ceos  (oi)+cos  (ir/4-a) 


and  therefore: 


k p |v  -V  |: 
w w 1 i w 1 


^wPw I vi I 

Cccos  (a)+cos  (ir/4-a)  J2 


we  have: 


tan 4*  » tan(rr/4-a)  + 


where : 


mfCCcoscH-cos  (TT/4-ot)  y 
kwpwco«  (Tr/4-a) 


Now  since : 


T = |*v|co«*  - Tix  - k^lv^rcoadr^-a) 
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we  get: 


k p | Vj  I 2cos  (TT/4-a) 
w w 1 


Lc cos  (a)  +cos  (ir/4-a)  U2 


, kwpw|  Vi  | 2cos  (it/ 4-a) 

T I . 

vi  

CCcos  (a)+cos  (ir/4-oO  ]2cos4‘ 


Now  since : 


we  get: 


kwpw^il2cos(7r/4'a) 

C^cos  (a)  +cos  (7T/4-a)  ]2cosHr 


or: 


|Vil 


|vj 


= C^cos  (a)+cos  (ir/4-a)  3 


kaPac°s  (HO 


kwPwcos  (TT/4-a) 


Now  from  above  we  have : 


J — = tan  (H1)  -tan  (tt/4-ci) 

|Vil 

and  multiplying  this  by  the  previous  equation  we  get: 


J? = Cscos  (ot)  +cos  (ir/4-g)  IL  / kaPaCos^ £tan  (¥) -tan  (tt/4-oi)  ] 

|va|  ^fkwpwcos  (tt/4-ci) 


The  equation  for  C is  derived  from  Ekman  theory,  as  follows: 


where:  Y = lv,_vwl 


y = the  coefficient  of  viscosity  between  ice  and  water 


i 


N|Vw|vf 


Where:  N = D/1v^|  and  is  the  depth  of  frictional  mixing.  This  has  been 

evaluated  as  N = 473,  (Shuleikin) . Now  since: 


C = — 


we  derive: 


V7Tk  V; 

Nf 


4 


7Tkwv2 


473f 


It  can  be  assumed  with  sufficient  accuracy  that  the  density  of  air  (pa)  is 
1.25x10  3 gm  per  cm3,  the  density  of  water  is  1 gm  per  cm3,  the  coefficient 
of  friction  between  ice  and  water  (kw)  is  4x10  3 and  between  air  and  ice 
(k  ) is  2xl0-3.  From  observations  it  has  been  determined  that  Ot  = 27°.  In 

3l 

addition,  since  the  density  of  ice  is  .9  and  the  mass  per  unit  area  of  ice 

(m)  is  the  density  times  the  thickness  (L) , m = .9L.  Therefore: 
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